Controlling and enhancing THz collective electron dynamics in superlattices 
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We show that a tilted magnetic field transforms the structure and THz dynamics of charge do- 
mains in a biased semiconductor super lattice. At critical field values, strong coupling between the 
Bloch and cyclotron motion of a miniband electron triggers chaotic delocalization of the electron 
orbits, causing strong resonant enhancement of their drift velocity This dramatically affects the 
collective electron behavior by inducing multiple propagating charge domains and GHz-THz current 
oscillations with frequencies ten times higher than with no tilted field. 
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Superlattices (SLs), comprising alternating layers of 
different semiconductor materials, provide a flexible en- 
vironment for studying quantum transport in periodic 
potentials and for generating, detecting, mixing, and am- 
plifying high-frequency electromagnetic radiation [THTT]. 
Due to the formation of energy "minibands", electrons 
perform THz Bloch oscillations when a sufficiently high 
electric field, F, is applied along the SL. The Bloch or- 
bits become more localized as F increases, thus produc- 
ing negative differential velocity (NDV) in the electron 
drift velocity, v^^ versus F characteristic |12l|13]. This 
single-particle NDV also strongly influences the collective 
behavior of the electrons, causing them to slow and ac- 
cumulate in high-density charge domains [T]-[3] . For suffi- 
ciently high F, the domains are unstable and propagate 
through the SL, generating current oscillations at GHz- 
THz frequencies determined by the form of Vd{F) and 
the SL length [HlliHlT]. In a given SL with fixed Vd{F), 
the frequency of domain dynamics can be tuned, over a 
limited range, by changing the applied voltage [6| fT4]. 

Here, we show that both the spatial profile of charge 
domains and their oscillation frequency can be flexibly 
controlled by using a tilted magnetic field, B, to engi- 
neer the Vd{F) characteristic of the SL. At F values for 
which the Bloch frequency equals the cyclotron frequency 
corresponding to the B component along the SL axis, the 
semiclassical electron motion changes abruptly from lo- 
calized stable trajectories to unbounded chaotic paths, 
which propagate rapidly through the SL p!8H22] . This 
delocalization creates a series of sharp resonant peaks 
in Vd{F)^ which were detected in previous DC current- 
voltage measurements p!9H2T1 , relate to mode coupling in 
Josephson junctions [23 , and can stabilize the SL Bloch 
gain profile in the vicinity of Stark-cyclotron resonances 
[24 . We show that these Vd peaks create multiple propa- 
gating charge domains, shaped by B and 6>, and thereby 
generate AC currents whose magnitude and frequencies 
are far higher than when B = 0. Chaos-assisted single- 
electron transport induced by the interplay between cy- 
clotron and Bloch motion therefore provides a mechanism 
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FIG. 1: (a) Semiclassical model in which the SL is described 
by a continuum region (dark gray) with miniband dispersion 
E(kx). We discretize this region into N layers. In the m*^ 
layer, the electron density is rim- At the left- [right-] hand 
edges of this layer, F = Fm [Fm+i]- Coordinate axes show 
orientation of F and B. (b) Vd versus r oc F calculated for 
B = 15 T with (from bottom to top) 6* = 0, 25° and 40°. 
For clarity, curves are offset vertically by 10^ ms~^. Inset: 
electron trajectories inx — z plane (scale bar at top) calculated 
over 4 ps for r = 1 (lower) and r = 7r/3 (upper) at ^ = 40°. 

for controlling the collective dynamics of the miniband 
electrons, thus increasing the power and frequency of the 
resulting current oscillations by an order of magnitude. 

We consider the GaAs/Al As/In As SL used in recent 
experiments [I9l|20]. Fourteen unit cells, each of width 
d = 8.3 nm, form the SL, which is enclosed by GaAs 
ohmic contacts [light gray in Fig. [ija)] with n-doping 
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density no = 10^^ m~^. Electrons are confined to the 
first miniband with l^inetic energy versus wave number, 
kx^ dispersion relation E(kx) = A[l — (zos{kxd)]/2^ where 
the miniband width A = 19.1 meV [19j. Semiclassical 
miniband transport corresponds to modeling the SL by a 
region of width L [dark gray in Fig. [ija)] where electrons 
move freely (with the GaAs effective mass m*) in the y—z 
plane but have dispersion, E{kx), along the SL axis. 

We calculated semiclassical trajectories for a miniband 
electron with B tilted at an angle to the SL (x) axis 
[Fig. [ija)] and F uniform throughout the SL [18l[19]. We 
then used an Esaki-Tsu (ET) model [12] [19] to determine 
= ('^a^(^) exp(— t/r))/ri, where Vx{t) is the electron 
speed along x at time t and (.) denotes averaging over the 
starting velocities of the hot miniband electrons, whose 
temperature 100 K) exceeds that of the lattice (4.2 
K), and integration over t, taking an electron scattering 
time r = ri[re/(re + ^i)]^^'^ = 250 fs determined from 
the elastic (interface roughness) scattering time Tg = 29 
fs and the inelastic (phonon) scattering time = 2.1 ps 
[201 [25] . When ^ = 0°, cyclotron motion in the y — z 
plane is separable from the Bloch motion along x. The 
lower curve in Fig. [TJb) shows Vd versus r = ujb/oJc ^ F 
where ujb = eFd/h is the Bloch frequency and ujc = 
eB cos is the cyclotron frequency corresponding to 

the X— component of B. As expected [lllHlllS], this trace 
peaks when r = I/ucT (i.e. ubt = 1) and thereafter 
decreases with increasing r as more electrons complete 
Bloch orbits before scattering. 

When ^ 7^ 0° , strong mixing of the cyclotron and Bloch 
motion drives the electron orbits chaotic [Fig. [TJb) in- 
set] [m [191 [21] . When r is irrational, the electron orbits 
remain localized along x [Fig. [TJb) upper inset]. By 
contrast, when r is an integer the electrons follow un- 
bounded paths [Fig. [ijb) lower inset] and map out in- 
tricate "stochastic web" patterns in phase space [181 US] • 
This abrupt delocalization of the electron paths generates 
sharp resonant peaks in Vd- The Vd{r) curves shown in 
Fig. [ijb) for = 25° (middle trace) and 40° (top trace) 
reveal a large peak at r = 1 and a smaller additional 
feature at r = 2, most apparent when = 40°. For r val- 
ues that are rational but not integer, the electron orbits 
are finite, but exhibit some resonant extension along x. 
This causes the small additional peaks visible at r = 0.5 
in the middle and top curves of Fig. [ijb). It is weh 
known that NDV in SLs and Gunn diodes creates prop- 
agating charge domains [1 . The multiple NDV regions 
associated with chaos-assisted resonant transport when 
^ 7^ 0°, suggests that the tilted B-field will induce more 
complex spatio-temporal domain dynamics. 

To investigate the collective behavior of the electrons, 
we solved the current-continuity and Poisson equations 
self-consistently throughout the device by adapting, for 
miniband transport, a model used previously to describe 
inter- well transitions in SLs [T] |l5l [IT]- In this model, 
we discretize the miniband transport region [dark gray 
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FIG. 2: I(V) characteristics calculated for (from bottom to 
top) ^ = 0°, 25°, and 40°. Curves are vertically offset by 
25 mA. Current oscillations occur within the shaded regions, 
whose upper [lower] bounds are -Trnax (V) [/min(V)]. Dashed 
curves are unstable steady state solutions of / 

in Fig. [ija)] into N = 480 layers, each of width 
Ax = L/N = 0.24 nm small enough to approximate a 
continuum. The volume electron density in the m^^ layer 
(with right-hand edge at x = mAx) is rim and the F val- 
ues at the left- and right-hand edges of this layer [vertical 
lines in Fig. [ija)] are F^, and Fm-\-i respectively. In the 
emitter and collector ohmic contacts, F = Fq. The evo- 
lution of the charge density in each layer is given by the 
current continuity equation 

eAx^^ = Jm-i - Jm^ m = l...N, (1) 

where e > is the electron charge and Jm = enmVd{Em)^ 
in which Fm is the mean field in the m^^ layer [27 , is 
the areal current density from the m^^ to the m + 1^^ 
layer neglecting diffusion [TTl [19] . Since Jm depends on 
the local drift velocity, Vd{Fm), the collective electron 
dynamics depend directly on the single electron orbits. 
In each layer, Fm obeys the discretized Poisson equation 

cAx 

Fm+i = {nm-nD)^Fm, m = l...N, (2) 

where Sq and Sr = 12.5 are, respectively, the absolute 
and relative permittivities and n^:) = 3 x 10^^ m~^ is the 
n-type doping density in the SL layers |19j . 

We use ohmic boundary conditions [l] to determine 
the current, Jo = aFo, in the heavily-doped emitter of 
electrical conductivity a = 3788 Sm~-^ [19]. The voltage, 
V ^ applied to the device is a global constraint given by 
V = U^^ EZ=i{Fm + i^m+i), where the voltage, /7, 
dropped across the contacts includes the effect of charge 
accumulation and depletion in the emitter and collector 
regions and a 17 1^ contact resistance [28]. We calculate 
the current I{t) = ^ T.l=o Jm^ where A = 5 x 10"^° 
m^ is the cross-sectional area of the SL [19 . 

Following initial transient behavior, I(t) either reaches 
a constant value or oscillates between minima and max- 
ima, /min and /max respectively, which depend on V, 5 
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FIG. 3: /(t) curves calculated for [0,V] = (a) [0°, 290 mV]; 
(b) [0°, 490 mV]; (c) [25°, 290 mV]; (d) [25°, 490 mV]; (e) 
[40°, 540 mV]; (f) [40°, 610 mV]. Arrowed peaks in (d) and 
(f) are discussed in text. Insets: Fourier power spectra, P(/), 
with a common vertical scale in arb. units and, for / > 0.2 
THz, also shown vertically enlarged xlO and offset. 



and 0. Figure [2] shows /(V) curves calculated for 5 = 15 
T and (from bottom to top) l9 = 0°, 25°, and 40°. Each 
trace is single valued at low but double valued in the 
shaded regions, where the lower [upper] boundaries show 
^min [^max]- For all ^, Stationary behavior occurs at low 
V where I{V) is approximately linear. But when V ex- 
ceeds a critical value, Vc^ which depends on B and ^, the 
stationary state loses its stability via Hopf bifurcation 
and I{t) starts to oscillate between /min and /max- The 
amplitude of the oscillations, la = /max — /min, increases 
with increasing V for all 0. In addition, for given V, 
la generally increases with increasing 0. In the regime 
where I{t) oscillates (shaded in Fig. |2|, there is also an 
unstable stationary state, found by setting drim/dt = 
in Eq. (1), corresponding to a fixed current whose V- 
dependence is shown by the dashed curves in Fig. |2j 
The shapes of the stationary I{V) curves, each compris- 
ing a stable (V < Vc) and unstable (V ^ Vc) part, are 
similar to previous experimental measurements \19\ . 
We now consider how the I{t) curves vary with V and 




FIG. 4: (Color) rim calculated for [0,V] - (a) [0°, 290 mV]; 
(b) [0°, 490 mV]; (c) [25°, 290 mV]; (d) [25°, 490 mV]; (e) 
[40°, 540 mV]; (f) [40°, 610 mV]. For clarity, upper (sur- 
face) plots are shown as grey-scale projections beneath where 
yellow, blue, purple and green curves are loci of constant F 
values corresponding, respectively, to the FT, r = 0.5, r = 1 
and r = 2 Vd{F) peaks. 



e. For l9 = 0° and V = 290 mV ^ K [Fig. [sja)], I{t) 
exhibits periodic oscillations whose frequency ~ 37 GHz 
corresponds to the single dominant peak in the Fourier 
power spectrum, P(/), inset. When V increases to 490 
mV [Fig. [sj^b)], the fundamental frequency of the oscil- 
lations falls to ~ 12 GHz. In addition, the peaks in /(t) 
sharpen, thus strengthening the higher frequency har- 
monics in P{f) [Fig. [Sjb) inset]. The I{t) and P{f) 
curves calculated for ^ = 25° and V = 290 mV ~ Vc 
[Fig. ^c)] are similar to those for l9 = 0° [Fig. ^d.)]. 
But when V increases to 490 mV [Fig. |3|d)], the results 
for = 25° differ markedly from those for ^ = 0° [Fig. 
[sjb)]. In particular, tilting B almost doubles the funda- 
mental frequency and introduces new I{t) peaks [arrowed 
m Fig. jsjd)] whose origin we explain below. Compared 
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with ^ = 0°, these extra peaks strongly enhance the high 
frequency components in P(/), as shown in Fig. [sj^d) 
inset, which reveals a dominant 3^^ harmonic at 54 GHz. 

When e = 40° and V = 540 mV ^ K [Fig. (sfe)], the 
shapes of the I{t) and P{f) curves are similar to those for 
^ = 0° and 25°, but the dominant fundamental frequency 
is higher. By contrast, when = 40° and V = 610 
mV [Fig. [sj^f)], the I{t) fluctuations are both stronger 
and richer than for l9 = 0° [Fig. (sfb)] and = 25° 
[Fig. [sjd)]. Consequently, the high frequency peaks in 
P(/) [Fig. [sj^f) inset] are greatly enhanced, with the 5^^ 
harmonic at 92 GHz being the strongest. 

To understand how I{t) varies with V and 6>, we con- 
sider how these parameters affect the underlying spatio- 
temporal electron dynamics. Figure [Ij^a) shows the sur- 
face plot of Umit^x) with its grey-scale projection be- 
neath calculated for l9 = 0° and V = 290 mV ^ Vc. For 
given X, rim oscillates periodically in t due to NDV in the 
corresponding Vd{r) curve [lower trace in Fig. 1(b)]. The 
dashed yellow curve in the lower projection of Fig. |4ja) 
shows the (t, x) locus along which F is fixed at the value 
corresponding to the ET peak in Fig. [ijb) lower inset. 
As X passes beyond this locus, the electrons slow, thus 
increasing the local values of both and F. This fur- 
ther decreases Vd and increases n^n, making the electrons 
accumulate in a charge domain [peak in surface plot in 
Fig. [4]^ a)], which propagates through the SL [1 . When 
the domain reaches the collector {x = L) it produces a 
sharp increase in I(t). Another charge domain then forms 
near the emitter and the propagation process repeats, so 
producing I{t) oscillations [1 . For larger V [Fig. [4]^b)] 
there are similar domain dynamics, but their frequency 
is lower because the higher mean F value reduces Vd- 

When V = 290 mV ^ increasing 6 from 0° [Fig. 
[4]^a)] to 25° [Fig. gc)] has little qualitative effect on 
the domain dynamics because V is low enough to ensure 
r < 0.7 through most of the SL: a regime where the Vd{r) 
curves for ^ = 0° and 25° have similar shapes [Fig. [ijb)]. 

This picture changes qualitatively when V becomes 
high enough to make r ^ 1 at some positions within 
the SL. Figure |4jd) illustrates this for V = 490 mV and 
= 25°. The yellow and purple curves in the lower pro- 
jection show the (t, x) loci along which F equals the val- 
ues corresponding, respectively, to the leftmost (ET) and 
r = 1 Vd peaks in Fig. [ijb). When t ^ 25 ps, NDV as- 
sociated with the ET peak creates a high density charge 
domain for x just beyond the yellow locus. At t ^ 50 
ps, a second charge accumulation region appears above 
the purple locus. This domain originates from the NDV 
region just beyond the r = 1 Vd peak. Its appearance 
produces an additional peak, labeled "r = 1" , in the I{t) 
trace in Fig. [sjd). When t ^ 65 ps, merger of the two 
charge domains creates the I{t) peak labeled "Merger" 
in Fig. [sj^d). After merger, the charge within the single 
domain is almost twice that for = 0°. In addition, the 
presence of the r = 1 Vd peak increases the mean electron 
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FIG. 5: (Color) Color map of fmax{V,0) for B = 15 T. 

drift velocity compared with 6> = 0°, thus also raising the 
domain propagation speed. These two factors increase 
both the frequency and amplitude of the I{t) oscillations 
[compare Figs. |3jb) and (d) and their insets]. 

Increasing to 40° further enriches the charge domain 
patterns. Since the r = 0.5, 1 and 2 resonances occur for 
smaller F at higher ^, their effect on the domain dynam- 
ics is apparent even for V very close to Vc. Figure [4]^e) 
reveals multiple charge domains near the yellow, blue and 
purple loci in the lower projection, along which F coin- 
cides, respectively, with the ET, r = 0.5 and r = 1 Vd 
peaks [upper curve in Fig. [ijb)]. Coexistence of multi- 
ple domains substantially increases both the amplitude 
and frequency of the I{t) oscillations [compare Figs, [sja) 
and (e)]. When V ^ 610 mV [Fig. [4]^f)] a new domain 
associated with the r = 2 resonance (green locus) ap- 
pears. The various domains produce multiple peaks in 
/(t), as shown in Fig. jsjf) where the labels mark peaks 
arising from formation of the r = 1 and 2 domains and 
their merger. These peaks create strong high-frequency 
components in P{f) [Fig. [sj^f) inset]. 

To quantify the effect of chaos-assisted transport on 
the I{t) oscillations, the color map in Fig. [s] shows the 
frequency, fmax^ corresponding to the largest peak in 
P{f)^ in the V — plane. For V < Vc (blue area in 
Fig. [5|, there are no charge domain oscillations. When 
V > Vc, where I{t) oscillations do occur, fmax generally 
increases with increasing 0, attaining a maximum (light 
yellow area) when V ^ 800 mV and ^ 70° . In this 
regime, fmax is ~ 10 times higher than for ^ = 0° due to 
the formation of multiple propagating charge domains. 

In conclusion, both the form and dynamics of travel- 
ing charge domains in a biased SL can be controlled and 
strongly enhanced by applying a tilted B-field. Addi- 
tional NDV regions created by B induce multiple charge 
domains, which increase both the amplitude and fre- 
quency of the oscillations in I(t) - giving an order of 
magnitude increase in the frequency and power of the 
dominant Fourier peak at large 0. Multiple Vd maxima 
can be created in other ways, for example by an AC elec- 
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trie field [TT] . Our results thus open routes to controlling 
the collective dynamics of charge domains in SLs by us- 
ing single- electron miniband transport to tailor Vd{F). 
It may also be possible to realize and exploit related dy- 
namics in nonlinear atomic and optical systems [29-31 . 
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